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The  Pile  Simulator: 


A  Network  Analogue  to  the  2-Group  Analysis.* 


1  The  Network  and  Its  Equations 

Consider  the  situation  shown  in  Fig.  1  in  which  a  current 

density  J  feeds  current  to  a  conducting  surface  (conductance  per  square 

of  at  a  point  whoso  potential  is  The  current  flows  away  (1)  in 

the  plane  with  surface  current  density  (2)  through  an  ardsotropic 

conductor  (•'specific'*  conductance  n^^  per  unit  area)  with  a  current  density 

to  a  second  surface  and  (3)  direct  to  ground  through  an  anisotropic 

condxictor  ("specific"  conductance  kn^^  per  unit  area)  with  a  current  density 

Ij.  The  current  is  dissipated  (1)  in  the  second  surface  with  current 

density  and  (2)  to  ground  through  an  anisotropic  conductor  ("specific" 

conductance  n2  per  unit  urea)  with  a  current  density  I2. 

The  equations  governing  the  electrical  relations  here  are: 

\  s  -  »  *2^^  ^2  =  '  *^2 

J  •  V  •  Kl  /  ll  /  I3  »  ^1  =  V*  ^  ^2 

‘ni  (?i-  Vg)  s  Ii  ,  Vj  s  I2 
ni  k  Vi  s  I3 


*  'Siis  analysis  has  been  made  in  coimection  with  the  design  of  a  pile 
slmuDAtor  using  point-by-point  ampUfioation.  The  first  such  device 
was  set  up  at  Oak  Ridge  by  E.  J.  ’Jtade  and  J.  Simpson. 


which  lead  to 


9  ,  “1  n-  1 

-  V  /  -  (1  /  k)  Vi  =  «2:  Vp  / _  J 


-V  ^2  / 


*2 


102 


®1 


(1) 


(2) 


k  neteork  analogue  to  Fig,  1  Is  shovm  in  Fig.  2.  The  meshes 
are  square,  S  ualte  on  a  side.  Accordingly,  the  current  fed  Into  euch 
upper  level  Jcnctlon  jnuet  be  JS^.  The  planar  network,  consisting  of  the 
conductances  g|^,  has  a  surface  resistance  per  square,  and  similarly 
for  the  lower  net.  Kccordlngly, 


«1  =  *1  »  82  -  “2 


The  interlevel  conductances  bj^  simulate  a  continuous  medliim  of  conductance 
2 

h^/S  ,  so  that 


"1  =  *  "2  = 


,2 


'ith  these  substitutions  Eqs.  (1)  and  (2)  become 


2  **1 

^  ''1  :s-  'i  •  -3^  '2  / 


^  81 


81 


& 


(3) 


.-2 


.  V  / 

2  -,2 


h^/h. 


Vo  : 


S‘  82 


dg2 


(4) 


II  The  Network  as  an  Knalogue 


In  the  pile  the  two-group  formulation  of  the  neutron  behavior 


-3 


-  D2  ^  2  ^2  ^2  =  ®2  >  intermediato 


T»here 


(f  =  ■>  »2,  , 


32:..r. 


^8l  ^Bl  ^ 


t2 


82 


Ao  s 


■s2 


B  -  ^ 

---..-  *  "2  *  "T" 


^4,5)>  (5) 


(6)  ,  (7) 


(8)  ,  (9) 


For  eaae  tn  satisfying  boundary  oondltiona  later  ><«  shall  use 
J  -  -J-  ^  5  ^8  ^  (defining  ;\g) 


30  that  the  group  equations  become 

2  T  .  ^3  J"  ^sl  ®1 


-  V  /  -i.  : 


^s2  ^1 


^82  ®2 


>'-/  f.  =irv^^ 

n2  D2 


(5) 


(6) 


Me  now  establish  the  correspondence  between  network  and  gro\^  equations, 
re  put* 

Vi  -  ®  >  ^2  ”  nyi/^2  (6e5) 


«  m  and  n  have  new  meanings  here  unrelated  to  ai  «  n^  ,  etc 


leaving  o  and ^  free  for  later  epociflcation.  In  the  simulator  JS^ 
idll  be  proportional  to  V2,  so  we  also  write 

nV2  an^ 

S2  S2 

»diere  n  ie  the  transfer  conductance  of  the  amplifier.. 

Bgr  equating  coefficients  it  follows  that 


(1  //(n  /  hjj  Bi 

81  h  '  32  8^^"  X2 


(10),  (11) 


S  g2 


D2 


“1  ^2  ®2 

/^S2g2  A]^  Dg 


(12),  (13) 


The  quantity  m  does  not  iq^ear.  Having  chosen  a  multiplying  median  as 
expressed  by  its  properties  in  the  ri^t  neabers  above,  having  set  k 
arbitrarily  for  the  time  being,  and  having  selected  the  scale  factor 
between  i^iper  and  lower  mash  voltages,  we  can  choose  a  conductance 
of  convenient  value  and  determine  in  turn  (from  10),  g2  from  (13), 
tram  (12)  and  n  from  (11)-  The  formal  analogy  between  network  and 
the  2>'grotqp  pile  equations  has  thus  been  established. 

Two  discrepancies  appear,  however,  idien  we  proceed  from  the  mathe^ 
matical  eynbols  to  the  physical  facts.  A  good  reflector  has  no  absorp¬ 
tion  and  no  fission.  This  means  that  63^  ^  A2  s  0.  The  first  dis¬ 
crepancy  lies  in  Eqo  (12)  which  requires  that 

h2  s  “  h^ 

which  is  inconvenient o  If,  actually,  we  make  h2  s  0,  then  the  network 
is  simulating  a  vli*tual  medium  in  which,  using  Eq.  (13) 


^.2 


s  Aa  = 


AZ  t 


The  virtual  Mditua  haa  an  absorption  rttleh  beooaos  mmQIst  if  /r  1b 
chosen  saallsr  and  ahleh,  naoerlcaUjr)  will  turn  out  to  be  nsgllfible,  bo 
that  this  virtual  nedluni  la  a  psmilsBlbls  sUbatltute  for  the  veal  nsdiUBo 
The  second  discrepancy  arises  frooBq,  (U>.  CoBd>lnlng  it  wltii 
Bq.  (10)  and  substituting  for  and  &]>  from  Sqs .  (6)  and  (7 )  leads  to 


n  ^f2  ^  ^  1 

Dl  /''^2  ^82  (1  i  k) 


(14.5) 


(1  /  k)  r,!  (ly  k)  ^ 


As  there  will  be  no  aiq>llfler  feed-back  in  the  reflector,  we  have 
n  s  0.  1st  us  achieve  this  in  our  virtual  reflector  by  letting 


»  Tn  =  0 


.1  ,,  ^ 

^,2 


This  virtual^  fission  is  again  a  property  which  is  under  control  thiou£^ 
Cortbinlng  Eqs,  (14)  and  (15)  gives 


1  /  k 


s2 


Sinc«  we  shall  find  that  k  is  of  the  order  of  unlty^  we  see  that  the 
virtual  reflector  yields  high  energy  neutrons  to  about  the  saos  extent 
that  it  absorbs  low  energy  oneso  'nius,  the  discrepancy  is  partially 
selfocoopensating  althou^  tending  to  Increase  the  fast  flux  and  to  de> 
crease  the  slow  flux  sli^tlyc  In  addition,  as  noted  above,  each  of  the 
virtual  properties  can  be  kept  small. 

Kith  the  preceding  discussion  in  mind  we  can  now  tabulate 
EqSo  (ID)  to  (13)  with  particular  reference  to  core  (subscript  C)  and  re¬ 
flector  (subscript  R)o  In  the  core  shall  have  no  need  of  direct  con¬ 
duction  from  the  high  network  to  ground  so  that  for  it  kg  -  0.  Ihen 


**1C  ^IC 

^SlD  °1G 


*^lli  ^  ^IR 

- - - -  -  -  (17J,  R) 

^  »ljl  ^ 


n  /  hiG  ^  Aic 

*3D  AA2C  Die 


D  -  0 


(18C,  R) 


hie  ^  B2e 

^  ®2C  ^IC  ^2C 


^IR  ^  AA2B  B2a 
s^g2R  Air  D2r 


(19C,  R) 


*^2C  ^  ^IC  t^2C 
°2C 


^2R  = 


(20G,  R) 


III  The  Boundary  Conditions  Between  Core  and  Reflector. 

The  boundary  conditions  between  core  and  reflector  which  are 
to  be  fulfilled  are: 


(21) 


»7 


(°(1,2)C  /^8(1,2)C^  ^^(1,2)C  =(°(l,2)R/\(l,2)R)'^i5l,2)B 


Between  core  net  and  reflector  net  the  conditions  that  are  autonatlealljr 
AilfUled  are: 

^(1,2)C 

^(1,2)C 

Equation  (21)  is  valid  by  virtue  of  Eq,  (6o5)(. 

Tm  eiiople  types  of  interface  between  networks  are  ^own  in 
plan  in  Fig»  3  at  A  and  Bo 

m  considering  the  second  boundary  condition  we  note  that  the 
electrical  analogue  to  grad  ^  is  grad  V,  idiich^  in  the  network,  becomes 
a  finite  voltage  difference «  In  order  that  it  apply  at  the  bouiklary, 
this  difference  should  be  taken  across  the  boundary  and  not  to  either  side 
of  it,  but  this  invades  the  medium  with  which  equality  is  to  be  establlshedo 
The  solution  of  this  difficulty  lies  in  Uie  stratagem  of  supposing  each 
region  independently  to  reach  across  the  boundary  just  for  the  purpose 
of  determining  the  change  in  V  vdilch  is  centered  at  the  boxmdary.  This 
is  Illustrated  in  Flgo  4A,  The  actual  boundary  arrangement  is  shown  in 
Figo  4Bc 

Considering  the  core  network,  it  is  seen  that 
(const)  X 


Hr 


2S 


K«  /  K" 

2d  sc 


(25) 


(1,2)B 


(1,2)R 


(23) 

(24) 


(22) 


By  splitting  the  vertical  conductors  into  halves  the  vertical  currents 


-e- 


are  halved  also  >riLth  the  residt  that  the  boundary  current  K  Is 

K«  K"  ^  ^ 

K  »  —11  n  I ..  —  -  (const)  X  wgQ 

2 

For  the  reflector  also, 

jj  -  (const)  X  3g|i  -  (26) 

^x 

and  idien  the  two  nets  are  cooblned  as  at  Flg«  4B  ee  have 


?  X  dx 


For  this  to  simulate  Eq*  (22)  m  nuist  have 

4(1,2)R 

«(l,2)R  ^(l,2)a  ^8(1,2)C 

This  establishes  the  relation  between  the  networks  representing  the  two 
OBdia  when  the  boundary  lies  as  at  Fig.  3Ac 

The  boundary  of  Pig.  3B  can  be  analyzed  In  the  sane  nay  and 
leads  to  the  sane  result. 

Figs.  3^'  and  3B'‘are  ailso  consistent  with  Eq.  (27),  the  only  dif¬ 
ference  In  the  analysis  being  that  coiqponents  of  the  flux  gradients  enter 
the  argument. 

The  explication  of  Eq.  (27)  to  Eqs.  (17C)  and  (17R)  establishes 
a  relation  between  the  conductance  levels  of  the  core  and  reflector  net- 


This  leaves  2  free  choices,  (1)  the  general  inpedance  level  of 
th®  networks  as  determined  by  the  initial  choice  of  gj^Q,  say,  and  (2)  the 
simulation  factor  ratio yV  .  The  first  is  deterjnined  by  the  requirements 
of  aiqsliflcation  and  sufficient  time  constant  and  tiie  second  by  the  need 
to  keep  the  virtual  absorption  in  the  core  (Eq,  ill*))  small <> 


IV  Interpretation  of  Critical  Condition  in  Network. 

Analysis  of  Eqa.  (17)  to  (20)  and  Eq.  (29)  reveals  that  a  value 
of  V  must  be  "buUt-into*'  the  network.  This  means  that  the  value  of  n 


frtiloh  sakes  the  siisalator  critical  cannot  be  interpreted  as  a  It 

can  be  interpreted  in  another  ivajTc  B^q  is  the  only  nuclear  property  co¬ 
efficient  which  involves  j  ^g2C  appears  only  in  Eq,  (18C) 

for  n.  We  can,  therefore,  look  upon  variation  of  n  as  effecting  vairiatlon 
of  the  low  velocity  fission  cross-section  in  accordance  with  Eq»  (l4o5) 


applied  to  the  core: 

Die  V/'^ac  Zsac 


(30) 


Another  approach  is  to  make  conprooises  ^ich  avoid  "building 
in"  V  o  This  requires  the  outright  assumption  that  s  0  which  makes 

independent  of  1) ^  This  would  mean  that  the  simulator  would  merely 
simolate  the  thermal  pile  case  vdth  adjusted  constants.  It  might  be  that 
relative  changes  in  would  have  fair  accuracy.  But  in  this  case  the 

V 

simulator  would  certainly  yield  a  bad  velocity  distribution  of  neutrons. 


7  A  Network  Design. 

For  an  octagonal  pile,  17"  between  parallel  faces,  with  18  unit 
msohee  in  that  distance,  S  s  s  2oA  cm.,  ^  5c76. 

For  the  nuclear  properties  R.  Ehrlich  and  1  have  used  the  Table 
of  Nuclear  Constants  of  CE-RE-1,  "A  Uultl-Grovp  Ifethod  for  Conputing 
Critical  Uasses  of  Intermediate  Piles,"  5/9/47.  CSroj;p8  A,  B,  and  C  were 
conbiaed  to  form  Gro\q>  1  here  and  0  and  E  gave  Group  2. 


-11- 


Ih«  tabulation  foUons: 


Core 

Beflector 

aroup 

Ai 

— 

D 

? 

or 

0 

2=. 

•5:. 

? 

■ 

H 

3.26 

3  27 

.0134 

.207 

.U5 

2.17 

1  67 

0 

207 

8 

2 

2cl6 

1.33 

.085 

0 

.180 

1.44 

.747 

B 

0 

(7) 

■ST,  r  (1 /=<)  2:t 

Tb«  0*«  of  Eqs.  (17)  to  (20)  appear  in  the  table«  For  the  other 

oonatanta  in  those  equations,  ne  then  have 

=  •0134  (1  -  )  /  —  =  -  .0157  /  o0259  5  o0102 

loU5  8 

Ajc  r  .085 

085 

s  2„5  —  s  .18  j  .0259 

XolB 

Ajjj  -  .0259  ;  ^2R  =  0  ;  r  0  j  z  .0259 
From  Eqc  (29)  k  =  1.54 

we  can  now  work  throu^  Eqs,  (17)  to  (20)  if  we  assune  an  hic  and  a  . 
let  h^Q  r  .25  x  lO”^  ,  -  4.0  jQ. ,  yUz  .1 

BlC  =  13.9  X  10“^  ,  =  72  K 

820  2  33.67  X  lo"^,  g2p“^  =  29.7  K 

h2Q  ;  11 .39  X  10  ,  h2(;  ~  87  o8  K 
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From  Eq.  (2?] 


g-  ”  10-66  X  10  ^  , 

XK 

«ia'^  “  93c9  k 

-1 

g2R  -  28  37  X  10  ^  , 

g2B  S  35-25  K 

hjjj  s  ,  376  X  ;  2  66 

(k  hjja)”^  ;  1.73 

Finally j  Eq  (14)  gives  for  the  virt''viJL  absorption  in  the  reflector 


_ aSS  r  ^00172 

^*e2R 

and  for  the  virtual  fission  Eq=  (16)  gives 


T" 

'-“82R 


.00068 


These  quantities  are  negligible. 


»»13~ 


VI  An  Atteopted  Generalization 

The  coireepondence  echeise  of  Eqo  (6o5)  ie  not  the  oioet  gmeral 
that  ni^t  be  used|,  and  it  haa  seemed  possible  that  by  generalizing  to 


z  au^l  “  »12  ^2 
$2  =  "®21^1  ®22  ^2 


(31) 

(32) 


the  need  for  keeping y<  small  in  order  to  minimize  ths  contradictions  of 
EqSe  (14)  and  (14«5)  mi^t  be  eliminatedo 

This  hope  has  not  been  fulfilled,.  BqSo  (31)  and  (32)  lead  to 
a  set  of  boundary  conditions  which  cannot  bo  satisfied  save  in  the  most 


exceptional  circumstances. 

Had  this  been  possible^  the  limitations  lii9)osed  b7  ^he  need  to 
have  ^  ^  X  v.ould  have  been  eliminated  and  it  might  have  been  possible 
to  cascade  several  groves  in  the  network.  As  it  lo,  the  practical  voltage 


limits  probably  make  it  iopractlcal  to  simulate  more  than  3  groiq>s<. 


\//^TOfiL _ 

£CT‘o^ 


T/^ua  _ . 

a  a  L  s.cT'o/^  ' 


Co/^6-  /ifi/O 
/?£at£CTO/? 


Figure  ka 

Analysis  of  Core-Reflector 
Boundary  and  Actual  Boundary  Arrangement 


